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We present results of calculations on Bose-Einstein condensed 87 Rb atoms subjected to a moving 
standing- wave light-potential of the form Vl(z, t) = Vo(t) cos(qz — uit). We calculate the mean-field 

C^j ' dynamics (the order paramter) of the condensate and determine the resulting condensate momen- 

tum in the z direction, P z (q, ui, Vo, t p ), where Vo is the peak optical potential strength and t p is the 
pulse duration. Although the local density approximation for the Bogoliubov excitation spectral 
distribution is a good approximation for very low optical intensities, long pulse duration and suf- 
ficiently large values of the wavevector q of the light-potential, for small q, short duration pulses, 

I '— ^ ■ or for not-so-low intensities, the local density perturbative description of the excitation spectrum 

breaks down badly, as shown by our results. 
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§ ' I. INTRODUCTION 

. 

S Elementary excitations of a Bose-Einstein condensate (BEC) can be explored with matter-wave interference studies 
using two-photon Bragg pulse spectroscopy [p]-pp[ . In such studies, the momentum imparted to the BEC by a Raman 
| scattering process can be studied as a function of the temporal duration of the optical pulses, t p , the detuning of the 
Bragg pulses from atomic resonance, A, the intensity of the Bragg pulses, the difference of the Bragg pulse wavevectors 
which is denoted by the wavevector q, and the difference of the central frequency of the two laser pulses ui. In the 
linear response limit, the response of the BEC to a weak perturbation with wavevector q and energy Tiw is given in 
terms of the dynamic structure factor 5(q, u>) flit], 
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C^' S(q,c) = I]Te-^" | (m\6p^\n) \ 2 8{uj - u mn ) , (1) 



£ — . where [3 — 1/kT, Tioj mn — E m — E n is the difference of two energy eigenvalues of the BEC, the density fluctuation 
Spq is induced by a perturbation with wavevector q and frequency u> that oscillates like e l (q r - wt ) an d z is the usual 
partition function. The momentum imparted to the BEC by the light-potential and its dependence on the wavevector 
q and frequency to can be directly related to the structure factor S^q, lu), and to the Bogoliubov dispersion relation, 
Eb (q) versus q (lj] . The excitation modes of a BEC have been measured Jl3|,[l4| and can be directly related to S^q, u>) 
and i?e(q). The momentum transferred to a BEC by a moving standing- wave light-potential has also been directly 
measured over a wide range of q and to [@,[|,[l0| , and it is of interest to calculate the momentum transfer versus q and 
'"O ' ^ so that the calculations can be compared with experiment. 

Here we report on the results of calculations of BEC excitation by Bragg pulse spectroscopic techniques to obtain 
the momentum transfer versus q and lo for weak and strong optical excitations. We find that, even with what is 
ordinarily considered weak intensity Bragg pulses, processes that are higher than first order (linear response) play 
a role in the excitation process. We also find that a local density Bogoliubov description is not valid for small 
wavevectors q. We describe the nature of the the higher order processes and their influence on the momentum and 
energy of the BEC excitation. We show how the simple Bogoliubov picture of the excitation is modified over a range 
of momentum-transfers and excitation energies due to higher order light-scattering processes, finite BEC size and 
inhomogeneity effects. 

We consider Bose-Einstein condensed 87 Rb atoms in the \F = 2, Mp = 2) hyperfine state confined in a harmonic 
oscillator potential, an array of optical traps and a gravitational field, and use parameters similar to those used in 
experiments carried out at the Weizmann Institute [ JL0| . The initial BEC is cigar shaped with N = 10 5 atoms in a static 

2 rnuj 2 — 

harmonic trap potential t4 Q (r) = ^^z 2 -\ ^ 2jL (x 2 +y 2 ) with frequencies lu z = 2ir x 25 Hz, u> x<y = uj x = uj y = 2ir x 220 

Hz (lu = (uOxUOyUOz) 1 ^ = 2tt x 106 Hz). The Bragg pulses propagate with wavevectors ki and k 2 in the x-z plane with 
angles ±9/2 relative to the x axis. The central frequency of one pulse is greater than that of the other, 0J2 = ^1 — 
and the difference frequency u is controlled using an acousto-optic modulator. The electric field takes the form 
E(£) = Ei(i) exp[z(ki -x— o7i<)+E 2 (i) exp[z(k 2 -x— uj 2 t), with lin || lin configuration for the field polarizations and equal 
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intensities for the two Bragg pulses. When the Bragg pulses are switched on, the atoms are trapped at the antinodes of 
a vertically oriented, red-detuned optical moving standing-wave; the antinodes are separated by Az — A/(2 sin(#/2)). 
The momentum transferred to an atom upon absorbing a photon from the field with wavevector ki and emitting a 
photon with wavevector k2 is given by ?iq = S,(ki — = hqz, where the Tiq = 2hk p h sin(0/2), and Tik p h = 2irh/\ is the 
photon momentum. Here A = 780 nm is the central wavelength of the Bragg pulses. The light-potential experienced 
by the atoms in the BEC as a result of the Bragg pulses is given by Vl (z, t) = Vo (f ) cos(qz — u)t) . The well-depth of the 
optical potential, Vo(t), is proportional to the intensity of the Bragg pulses and inversely proportional to the detuning 
from resonance, A, i.e., Vb(i) = MhJJghMl where Oi (£) and ^(f) are the Rabi frequencies. The well-depth temporal 
dependence /(f), where Vb(f) = Vo /(f), is taken to have a Gaussian rise-time and fall-time of width t r = 20 fj,s; Vb(f) 
is constant for a time duration f p between the rise and fall (/(f) = 1 for the time interval t p so Vb(f) = Vo in this 
interval). Pulses with short (t p = 1 ms) and long pulse duration (t p = 6 ms) are used. The strength of Vo that is 
used in the calculations will be specified in recoil units, En = (hkph) 2 / (2m) . Absorption of a photon from one pulse 
and stimulated emission into the other pulse produces a perturbation with energy fiw and momentum Tiqz. In the 
experiments, the light-potential and the harmonic potential are both switched-off (dropped), releasing the atoms to 
fall under the influence of gravity, and absorption images are taken after the particles evolve for some specified period 
of time under the influence of gravity. 



II. THEORETICAL FORMULATION 



We calculate the dynamics within a mean-field treatment using the time-dependent Gross-Pitaevskii equation (GPE) 

at 2m 

where 

V(r,t) = V h o(r)-mgz + V L (z,t), (3) 

is the potential, and g — 47TN £° h is the atom-atom nonlinear interaction strength which is proportional to the s-wave 
scattering length ao and the total number of condensed atoms N. The wave function (order parameter) ip(r, f) is 
propagated with a split-operator fast Fourier transform method. Due to the large number of grid points necessary 
in the lattice direction (z), we found it necessary to convert the 3D GPE into an effective ID GPE with similar 
dynamics as described in Ref. Jl7| . The wave function in momentum space, ip(k z ,t), is determined by taking the 
Fourier transform of ip(z, t). The net momentum of the BEC at any time f is given by 

/oo 
dk z r(k z ,t){hk z )4>(k z ,t). (4) 
-oo 

The expectation values of all dynamical quantities (e.g., energies, (Az), etc.) can be easily determined using the 
calculated wavepackets in either position or momentum space. 

A perturbative estimate of the rate of momentum-transfer to the BEC by the Bragg pulses is given by Jul 

dP z {t) 2n fV \ 2 

-^r = Y q {Y) [S(q^)-S(-q,-^)]. (5) 
For a uniform zero-temperature BEC, the Bogoliubov excitation energy at momentum-transfer q is 



Hence, the momentum transferred is given by 



E B (q) = hn B (q) = y/e(qMq) + 2gn] , (6) 

e(g) , 
E B (q) 1 



where e(q) — (hq) 2 /2m, n is the density, and the dynamic structure factor is given by 5(q, uj) — ; 5(hu> — Es{q)) 



Thus, the momentum transferred by a moving standing-wave excitation to a uniform condensate is non-vanishing 

Eb(<j)- 



only when u> = fig (q) , and its magnitude is proportional to the product of the integral over time of V Q 2 and q< ' " ' 
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For a non-uniform BEC with density profile that varies smoothly with position, one can define a local Bogoliubov 
excitation energy, 



E B (q,r) = hQ B (q,r) = ^e{q)[e{q) + 2gn(r)} , (8) 
and a local density approximation dynamic structure factor that behaves locally as a uniform gas [Q: 

Slda(ci,lu) = N- 1 J drn(r) E ^ 5(fiuj - E B (q, r)) . (9) 

Hence, in perturbation theory (i.e., small Vo), the momentum transferred to the non- uniform condensate can be 
approximated by substituting Slda(ci, w) for S(q,u>) on the right hand side of Eq. (|^). Thus, momentum transfer 
of a non-uniform condensate via a moving standing-wave excitation is smeared over a range of frequencies around 
10 = fl B (q), and its magnitude is determined by carrying out an average of E q ^\^ over the local density. 

III. RESULTS 

Calculations were carried out over a range of momentum-transfer q (i.e., a range of angles &), frequencies uj, pulse 
duration times t p , and laser pulse intensities (potential strengths Vo). 

Fig. [I] shows the calculated momentum imparted to the BEC by the Bragg pulses, P z (q, £*;, Vo, t p ), as a function of 
w for 6 = 22.5° (q = 0.390 %k p h — 80\/hmuj z ) and pulse duration t p = 6 ms. The calculations were performed for 
three different laser intensities corresponding to potential well depths Vo = 0.0054, 0.0162 and 0.054 Er. The low 
intensity curve peaks near the Bogoliubov frequency Cl B (?) = 8657 s _1 . This peak has a tail at lower uj which is due 
to the non-uniform density of the BEC; E B (q,r) for positions away from the center of the BEC is smaller than at 
the center and this can give rise to the tail, as can be understood from Eqs. (f|) and (||). At intermediate and high 
Vo, a power broadening of the spectral distribution P z (uj) is evident in Fig. Of; at these values of Vo, higher order 
(nonlinear) processes that populate ±2g momenta take place (see Fig. ^|). These phenomena can not be understood 
from a perturbative treatment. At even higher values of Vq, the spectral distribution becomes even wider and the 
peak structure becomes even more complicated and ragged. 

Fig. H shows \ip(k z ,t)\ 2 vs. k z for the low intensity case appearing in Fig.[l] for uj = 8945 and 9597 s _1 . The 
net positive momentum, (P z (uj) > 0) resulting at these frequencies in Fig. [yis due to the fact that the peak at 
k z = q (w 0.390 hk p h = 80\/hmuj z ) is larger than that at k z = —q. The peak near k z — q has an additional high- 
frequency feature at around k z — l20\/hmuj z ; the origin of this feature is not clearly understood. The peaks at 
k z = ±2q are almost two orders of magnitude reduced compared with the k z — ±q peaks. As the intensity of the 
lasers (the potential well depths) increase, the size of the k z — ±2q peaks grow in comparison with the k z — ±q peaks, 
and \ip(k z ,t)\ 2 grows at intermediate values of k z between the peaks at k z — ±jq, where j is a positive integer. 

Fig. |^ is similar to Fig. [j], except the pulse duration was taken to be 1 ms. The width of the distribution P z (uj) as 
a function of uj for t p = 1 ms is considerably wider than for t p = 6 ms for the low and intermediate intensities. The 
shorter temporal duration and therefore larger bandwidth of the 1 ms Bragg pulses allows for a wider distribution 
in frequency of P z (to) versus uj. The width hides the tail of the distribution at lower uj due to the non-uniform BEC 
density. For the high intensity case, the distribution does not change very much from the low and intermediate cases. 
Moreover, it is narrower than the high intensity 6 ms result shown in Fig. |l| because the pulse fluence is smaller, hence 
higher order processes do not significantly broaden the distribution. However, for the ultrahigh intensity case, power 
broadening of the distribution is significant. Note that the perturbation theory expression (||) can not account for the 
time-domain broadening shown in Fig. ^ (since it is derived assuming that the spectral width of Vq is within that of 
S(q, uj)). Eq. (||) shoud be modified to account for the bandwidth of the optical pulse: 

dPz[t) ~ f4= / *|VoM| 2 S(q,a,) + c.c) . (10) 



dt 2Ti VV2^ 

This equation should then be integrated over time to obtain an expression that replaces (^) for P z (t) . 

We now consider smaller values of momentum-transfer q such that the excitation wavelengths of the optical potential 
are comparable or larger than the size of the condensate. Fig. || shows P z versus uj for 9 = 10° (q = 0.174 Tik^ = 
S5.8^/hmuj z ) and t p = 6 ms. The curve for low intensity, Vo = 0.0054 Er, is similar to the low intensity result in Fig. ^ 
in the sense that a peak exists near the Bogoliubov frequency Cl B (q) — 4164 s _1 and a low frequency tail is present. 
Increasing the intensity to Vo = 0.054 Er (the curved labeled high), yields a multi-peaked ragged spectrum due to 
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higher order nonlinear processes. Upon reducing the optical potential by a factor of 0.3 from the low intensity case 
to Vq = 0.00162 En (the curve labeled ultralow), new features in the spectrum become clear. The small feature near 
lu = 860 s _1 and the tail of the Bogoliubov peak at 2411 s _1 , which was presumably due to an inhomogeneous density 
effect, become two ancillary peaks. The reason for the structure in the spectrum for the ultralow (and low) intensity 
can be understood by looking at \ip(k z ,t)\ 2 vs. k z for this case as shown in Fig. ||. Peaks at k z = q = 35.8\/hmuj z 
are clearly seen for ui — 3807 s _1 and these peaks fall within the tail of the oscillatory structures associated with the 
broadened k z = condensate. These structures at k z = q are responsible for the peaks in the spectral distribution 
in Fig ^. The peak near uj — 860 s _1 and the minimum near lu — 1636 s _1 in Fig. ^ result due to subtle interference 
of the k z — q peak with the structure surrounding the central (k z — 0) peak in Fig. || This kind of interference can 
occur when q (angle 9) is sufficiently small that the peak at k z — q is within the structure of the central peak. It does 
not occur at 9 = 22.5° (q = 0.390 hk p h), rather only for angles 9 < 10°. 

In Fig. ^ we plot P z versus u> for 9 = 5° (q — 8.72 x 10~ 2 hk p h = 18 y/hmcu z ), t p — 6 ms and potential well depths 
Vq = 0.00162, 0.0054, 0.0162 and 0.054 Er. At ultralow, low and intermediate intensities, an additional peak appears 
at lu w 400 s^ 1 . At all but ultralow intensity, the main peak in the spectrum is near the Bogoliubov frequency 
Qb(q) = 2061 s _1 . At ultralow intensity the peak at lu « 400 s _1 is even larger than that at lu s» 2000 s _1 and 
P z between these peaks becomes negative. Again, interference effects arising for reasons explained in connection 
with Fig. are apparently responsible. A complicated interference pattern appears around the values of k z = ±q in 
\ip{k z ,t)\ 2 vs. k z (figure not shown). This interference plays a role in the determination of the spectrum shown in 
Fig. The peaks in |^(fc 2 )| 2 at k z = q for w — 370 and 1687 s _1 are larger than the peak at k z = q for lu = 990 s _1 , 
and therefore maxima occur in the distribution shown in Fig. ^| at lu — 370 and 1687 s^ 1 and a minimum between 
the peaks in the spectrum occurs for lu — 990 s _1 . 

Fig. shows P z versus lu for 9 = 3° (q = 5.24 x 10~ 2 %k ph = lly/hmuj z ), t p = 6ms and V = 0.00162, 0.0054, 0.0162 
and 0.054 Er. The peak near lu — Hb(q) — 1234 s _1 for intermediate and high intensities becomes broadened at low 
intensity and then splits into two peaks at lu w 400 and 1500 s _1 for ultralow intensity, with the one at lu w 400 s _1 is 
about four times larger in magnitude than the one atw« 1500 s _1 . The dip between the two peaks for the ultralow 
intensity is apparently also an interference effect that can not be explained in terms of a local density approximation 
to the Bogoliubov spectrum. For intermediate and high intensities, the peaks in \i/j(k z )\ 2 at k z — q are sufficiently 
large that interference with the structure around the k z = peak does not occur and therefore the interference dip is 
absent. 

IV. SUMMARY AND CONCLUSIONS 

For very low optical intensities, long pulse duration and sufficiently large values of the momentum-transfer imparted 
by the light-potential, the local density approximation for the Bogoliubov excitation spectrum is a reasonable first ap- 
proximation; the response peaks in P z (g, u>, Vq, t p ) versus lu are broadened, particularly to lower values of lu. However, 
we have shown that even for relatively low optical intensities, power broadening results and higher order processes 
occur that correspond to moving the atoms from the wave packet with central momentum k z = q to wavepackets 
with central momentum 2q and as well as additional k z = nq, with n > 2 and n < —I. Moreover, at lower values 
of momentum-transfer q (smaller angles 9) where the wavelength of the optical potential becomes comparable to or 
larger than the size of the condensate, interference effects play a role in the dynamics, and directly affect the spectrum, 
P z (q, lu, Vq, t p ) versus lu and can produce additional maxima and minima in the spectrum. The power spectrum of the 
order parameter after the optical potential is dropped, \ip(k Zl t p )\ 2 versus k z , can be used to understand the nature of 
the spectra P z (q,tu, Vo,t p ) versus lu and q. 

One should keep in mind that the study performed here is a ID calculation of the full 3D dynamics; 3D effects may 
modify details of the results we obtained. Nevertheless, we believe that the qualitative features of the conclusions 
will not change. We have detailed elsewhere how our quasi- ID calculations of the type we presented here model 3D 
aspects of the dynamics in cylindrically symmetric potentials pl[ , but this method can not describe radial excitations 
of the BEC that might arise due to the optical potential via the mean-field interaction. To the extent that radial 
excitations are not important, our method should be an adequate approximation to the 3D dynamics. 
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FIG. 1. Momentum P z (q, ui, Vo, t p ) versus uj for 6 = 22.5° (q = 0.390hk p h) and pulse duration t v = 6 ms. Results for three 
different laser intensities corresponding to potential well depths Vo = 0.0054 (low), 0.0162 (intermediate) and 0.054 (high) Er 
are shown. 
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FIG. 2. j?/>(fc z ,t)j 2 vs. k z for the low intensity case appearing in Fig. |l] at 
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FIG. 3. Same as Fig. [jj except the pulse duration is 1 ms and an additional well depth of Vo = 0.54 (ultrahigh) Er is 
included. The arrow indicates the position of f2_g(g). 
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